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1. Introduction 

An observation of the lepton-flavor-violating (LFV) or the lepton-number-violating (LNV) 
processes is a clear signature of the existence for physics beyond the standard model (BSM) since 
they are forbidden in the standard model (SM). Among various LFV processes, some LFV decays 
of tau lepton are expected to be enhanced by several kinds of BSM, especially, models with su- 
persymmetry. One of the most frequently discussed is x~ — > jx^y decay. But, in some models, 
such as the non-universal Higgs mass model, the constrained minimal supersymmetric SM and so 
on, among the X LFV decays, x~ — s> jJL tj or x~ — s> /i~p° will be enhanced. [|lj] In this talk, we 
report the recent result of a search for T~ -> £~M° (£ = e,H, M° = n°,ri,ri',p ,K* ,K* , co and <p) 
with the world largest data sample accumulated at Belle. On the other hand, some LNV decays are 
possible when the neutrino is the Majorana type spinor. Here, we show the first result of search for 
B+ -> D-£+£+ at Belle. 

2. t — > £M° {£ = e,n, M° = ^°,T],r] / ,p ,K*°,K* , (o and <j>) 

2.1 Analysis Method 

In the x LFV analysis, in order to evaluate the number of signal events, two independent 
variables are defined, that are signal-reconstructed mass and energy in the center-of-mass (CM) 
frame from energies and momenta for the signal x daughters. In the x — > jUTj case, they are defined 
as = ^Je~^~-P~^, AE = E^ — E^ a , where (P^) is a sum of the energies (a magnitude 
of a vector sum of the momenta) for pL and tj, the superscript CM indicates that the variable is 
defined in the CM frame and i^eam means tne initial beam energy in the CM frame. Principally, 

and AE should be m x (~ 1.78 GeV/c 2 ) and (GeV), respectively, for signal events while 
and AE will smoothly vary without any special structure in the background (BG) events. Due to a 
finite resolution, the signal events are distributed around ~ ra T and AE ~ (GeV). Taking into 
account the resolution, we set the elliptic signal region. Finally, we evaluate the number of signal 
events in the signal region. When the number of the observed events is consistent to that of the 
expected BG events, the upper limit for the number of the signal events is evaluated by Feldman- 
Cousins method. |Q] To avoid any bias, we perform the blind analysis: Before fixing the selection 
criteria and the evaluation for the systematic uncertainties, we cover the data events in the signal 
region. 

2.2 x^£P° (£ = e,n,P Q = 71°, 17,77') 

We perform a new search for the x decay into a lepton (e or }x) and a neutral pseudoscalar 
(71°, T7 or r]') with a 901 fb _1 data sample. A neutral pion is reconstructed from 2 photons while 
an tj (tj') is reconstructed from 77 (p°7) as well as (7T + 7r~Tj) to increase the detection 

efficiency. When the neutral pseudoscalars are reconstructed, their four-momentum is evaluated by 
a mass-constrained fit to obtain a better resolution for the signal region. Because we have modified 
the selection criteria applied to the previous analysis, we obtain an about 1.5 times better detection 
efficiency while the similar background level is kept. As a result, we observe one event in the x — > 
eTj (— > 77) mode while no events are found in other modes as shown in Fig. [jj Since these results 
are consistent with the background estimation, we set upper limits on the following branching 
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Figure 1: Resulting 2D plots for x — >• eft (a), e?](— s- 77) (b), er]{-^ nnn ) (c), eT}'(— > 7P°(— > ftft)) (d), 
7Z7zr](-> 77)) (e), ^7T° (f), jJ.r](-^ 77) (g), ^rj(^ 7r7r7r°) (h), A"?'(-> 7P°(-> tttt)) (i) and £trj'(->- 
7C7rr}(— !> 77)) (j), on the M£po - A£ plane. Here, black dots (shaded boxes) express the data (signal MC), 
the region bounded by two lines is defined as a 3(7 band for the BG estimation and the elliptic region is the 
signal one which corresponds to 3(7 in each plot. One event is found in the signal region for x — > et] while 
no events appear in any other modes. 



fractions: ^(t~ -> e~n°) < 2.7 x 1(T 8 , £ 

i<r 8 , #(t~ -> ii-ti) < 2.3 x ifr 8 , m{%- 

10" 8 , at the 90% confidence level. 



(%- -> jj.-K°) < 2.2 x ifr 8 , 

-> e _ Tj') < 3 - 6 x 10_8 and 



-> <? 77) < 4.4 x 
A* r? 



') < 3.8 x 



2.3 t -> ^° = e, V° = p°,K*°,K*°, co, 0) 

Similarly to t~ — > £~P°, we update our results of the search for t~ — > £~V°, where £ = e,{X, 
V° = p°,K*°,K* , ou, <p with an 854 fb _1 data sample. By performing a detailed background study, 
we obtain a 1 .2 times better efficiency in average with keeping similar level backgrounds. Finally, 
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Figure 2: Resulting 2D plots for x ep° (a), eK*° (b), e£*° (c), ecu (d), eQ (e), Atp° (f), jlK*° (g), At# *° (h), 
|i(0 (i) and ji<p (j) on the M f:v o - AE plane. Here, black dots (shaded boxes) express the data (signal MC), 
the region bounded by two lines is defined as a 5(7 band for the BG estimation and the elliptic region is the 
signal one which corresponds to 3<7 in each plot. One event is found in the signal region for x — > jJ.K*°, 
fiK*° and fj,<f> while no events appear in any other modes. 
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one event is found in the signal region for T~ — > n~K*°, jj,~K*° and while no events appear in 
any other modes. They are consistent with the expected number of the backgrounds. Consequently, 
we set the 90% confidence level upper limits on the branching fractions: £3{x~ — > e~p ) < 1.8 x 
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3. B+^i+H+D- (£,£' = e,n) 

If the neutrino is a Majorana-type spinor, a 8 + 4 
£ + £' + h~ decay is possible, where £,£' = e,}l, h = 
7t, K,p,K*,D, ■ ■ • [g]. Due to the size of CKM matrix, 
B + — > l + £' + D~ is expected to have the largest branch- 
ing fraction among them. But this mode has never been 
measured yet. Using a 7.7 x 10 8 B + B data sample, we 
perform a first search. Similarly to the X LFV case, we 
define M^ c and AE, where they are reconstructed from 
£ + , £' + and D , but the total energy is set to the initial 
beam energy. In order to evaluate the number of signal 
events, the signal box is defined as 5.27 (GeV/c 2 ) < M\, c 
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Figure 3: Decay process for B + — > 
£ + £' + D~ . When a neutrino is a Majorana 
spinor, the neutrinos can connect as an in- 
ner line in the process. 



< 5.29 (GeV/c 2 ) and -0.035 (GeV) < AE < 0.035 (GeV) for B + -> ii+ij+D , or -0.055 (GeV) 

< AE < 0.035 (GeV) for others, because the distribution for the electron energy has a small tail 
on the lower side. After the selection, we found no events in the signal box for all modes. Finally, 
we set the 90% confidence level upper limits on the branching fractions: £$(B + — )■ e + e + D~) < 
2.7 x 10~ 6 , S8(B + -> e + ix + D-) < 1.9 x 10~ 6 , m{B + ii+ii+Cr) < 1.1 x 10~ 6 . 
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Figure 4: Resulting 2D plots for B+ -4 e+e + D~ (a), B + -4 e+ji+D- (b) and B + -> fi+jj, + D- (c) on the 
- AE plane. Here, black dots express the data and the box region is the signal one in each plot. No 
events are found in each mode. 
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